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ABSTRACT 

We present the first images of the J=5-4 and J=16-15 lines of the dense gas tracer, cyanoacetylene, 
HC3N, in an external galaxy. The central 200 pc of the nearby star-forming spiral galaxy, IC 342, 
was mapped using the VLA and the Plateau de Bure Interferometer. HC3N(5-4) line emission is 
found across the nuclear mini-spiral, but is very weak towards the starburst site, the location of 
the strongest mid-IR and radio emission. The J=16-15 and 10-9 lines are also faint near the large 
HII region complex, but are brighter relative to the 5-4 line, consistent with higher excitation. The 
brightest IIC3N emission is located in the northern arm of the nuclear minispiral, 100 pc away from 
the radio/IR source to the southwest of the nucleus. This location appears less affected by ultraviolet 
radiation, and may represent a more embedded, earlier stage of star formation. HC3N excitation 
temperatures are consistent with those determined from C^^O; the gas is dense, 10'*"^ cm~^, and 
cool, Tfe < 40 K. So as to not violate limits on the total H2 mass determined from C^^O, at least 
two dense components are required to model IC 342's giant molecular clouds. These observations 
suggest that IIC3N(5~4) is an excellent probe of the dense, quiescent gas in galaxies. The high 
excitation combined with faint emission towards the dense molecular gas at the starburst indicates 
that it currently lacks large masses of very dense gas. We propose a scenario where the starburst is 
being caught in the act of dispersing or destroying its dense gas in the presence of the large HII region. 
This explains the high star formation efhciency seen in the dense component. The little remaining 
dense gas appears to be in pressure equilibrium with the starburst HII region. 
Subject headings: galaxies: individual (IC 342) — galaxies:starburst — galaxies: ISM — radio lines 



1. INTRODUCTION 



Little is known about the properties of the dense gas 
component on CMC scales in external galaxies yet it is 
dense gas (n//^ > 10^ cm"'^) that most directly correlates 



with global star formation rates (e.g., [Gao fc Solomon 
2004||Wu et~al]|2005[ ). To understand the evolution and 
regulation of current star formation in galaxies, dense 
gas properties need to be accurately characterized. 

In a recent imaging survey of millimeter-wave molecu- 
lar lines in the nearby spiral nucleus of IC 342, a surpris- 
ing degree of morphologica,! vari ation is observed in its 
dense gas ( Meier fc Turner|[2005), which was not antici 



pated from CO or HCN (e.g. , Ismzuki et al. 1990 Downes 
et al. I992 I). These and other observations (e .g 



Garcia- 



gurilio et al]l2000| [20011 12002| [Usero et al.l[2004, 2006^ 
conhrm that CO and MCN fail to fully constram excita- 
tion and chemical properties of the dense clouds. Various 
probes of the dense component exist, each with various 
strengths and weaknesses (e.g., Meier & Turner '2005'; 



Papadopoulos|2007, .Bussmann et al., 2008, .Narayanan et 


al.| 


2008 Krips et al. 2008 Gracia-Carpic 


et al. 


2008). 



primarily because it is bright, however, it is optically 
thick in most st arburst envi ronments (e.g . Downes et al. 
19921 iMeier & Turner 2004j iKnudsen et al. 2007 ), has a 
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widely spaced rotational ladder with its J>3 transitions 
in the submm and can be sensitive to chemical effects 
(PDR effects and IR pumping; e.g 
Aalto et al.||2002 ). 



Fuente et al. 1993 



Mere we present a study of the physical conditions 
of the very dense gas component of a nearby starburst 
nucleus using HC3N. With its large electric dipole mo- 
ment (^ = 3.72 debye versus 3.0 debye for HCN), low 
opacity and closely spaced rotational ladder accessible 
to powerful centimeter and millimeter-wave interferom- 
eters, HC3N is well-suited to high resolution imaging of 
the structure and excit ation of the densest component of 



et al. 1983 



the niolecular gas (eg. Morris et al. 1976 Vanden Bout 



2OIIL HC3N(5-4) has a crit- 
■3 — - 



Aladro et al 

ical density of 5 x 10* cm~'^ and an upper level energy 



of E/k = 6.6 K; the 16-15 line has a critical density of 
5 X lO'^ cm~'^ and an upper level energy of E/k = 59 K. 

We have observed the J = 5-4 and J — 16-15 lines of 
HC3N at <2" resolution in IC 342 with the Very Large 
ArrajF] (VLA) and the Plateau de Bure Interferometer 
(PdB)7 respectively. IC 3 42 is one of the c l osest (D ^3 



sev 



Mpc, or 2" = 29 pc; eg. |Saha et al.]|2002| [Karachent 



|2005|), large spirals with active nuclear star forma- 



tion ( [Becklin et al.|1980| [Turner fc Ho|1983^ . It is nearly 
face-on, with bright molecular li ne emission both in dense 
clouds and a diffuse medium ( L o et al. '1984; Ishizuki[ 



et al.[[1990[ [Downes et aL][1992f[Turner fc: Hurt 1992^ " 



I'he wealth of data and proximity of IC 342 permit the 
connection of cloud properties with star formation at 
sub-CMC spatial scales. The data pubhshed here rep- 

^ The National Radio Astronomy Observatory is a facility of 
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resent the first high resolution maps of HC3N(5-4) and 
HC3N(16-15) in an external galaxy. Excitation is a key 
component in the interpretation of molecular line inten- 
sities. So these two maps are compared with the previ- 
ously published lower resolution HC.-jNf lO-Q) map made 
with the Owens Valley Millimeter Array ( [Meier &: Turner 
2005 1 , to constrain the physical conditions of the dens 



est component of the ISM in the center of IC342, and to 
correlate properties of the dense gas with star formation, 
diffuse gas, and chemistry. 

2. OBSERVATIONS 

Aperture synthesis observations of the HC3N 3=5-4 
rotational line (45.490316 GHz) towards IC 342 were 
made with the D configuration of the VLA on 2005 
November 25 (VLA ID: AM839). The synthesized beam 
is l'.'95xl'.'5 (FWHM); position angle (pa)=-11.2°. Fif- 
teen 1.5625 MHz channels were used, for a velocity res- 
olution of 10.3 km s~^, centered at vlsr = 30 km s~^. 
The phase center is a(J2000) = 03'^46'"48!3; (S(J2000) = 
68°05'47'.'0. Amplitude tracking and pointing was done 
by observing the quasar 0228-1-673 every 45 minutes. 
Phases were tracked by fast switching between the source 
and the quasar 0304-1-655 with 130s/50s cycles. Absolute 
flux calibration was done using 3C48 and 3C147 and 
is good to ~5-10%. Calibration and analysis was per- 
formed with the NRAO Astronomical Image Processing 
Software (AIPS) package. The naturally weighted and 
CLEANed datacube has an rms of 0.70 mJy beam~^. 
Correction for the primary beam ('-^l' at 45 GHz) has 
not been applied. The shortest baselines in the dataset 
are ~5 kA, corresponding to scales of ^40"; structures 
larger than this are not well-sampled. 

The HC3N(16-15) line emission in the center of IC 
342 was observed with 5 antennas using the new 2 mm 
receivers of the IRAM Plateau de Bure interferome- 
ter (PdBI) on 2007 December 28 and 31 in C config- 
uration with baselines ranging from 24 to 176 m. The 
phase center of the observations was set to a(J2000) = 
03'>46'"48'rl05; (S(J2000) = 68°05'47'.'84. NRAO 150 and 
0212-1-735 served as phase calibrators and were observed 
every 20 min. Flux calibrators were 3C84 and 3C454.3. 
The calibration was done in the GILDAS package fol- 
lowing standard procedures. The HC3N(16-15) line at 
145.560951 GHz was observed assuming a systemic ve- 
locity of vlsb. — 46 km s^^ and a spectral resolution of 
2.5 MHz (5.15 km s~^). The average 2 mm continuum 
was obtained by averaging line-free channels blue and 
redward of the HC3N and H2CO lines and subtracted 
from the uv datacube to obtain a continuum-free dat- 
acube. The final naturally weighted and CLEANed dat- 
acube with lOkms"^ wide channels has a CLEAN beam 
of 1.83" xl.55", pa=46° and an rms of 1.7 mJy/beam. 

Single-dish observations of HC3N(16-15) find a peak 
brightness temp erature of 6.7 mK in a 16.9 beam (Al- 
adro et ar||2011 ). Convolving our map to this beam and 
sampling it at the same location yields a peak bright- 
ness of 7.1 mK, which agrees within the uncertainties of 
both datasets. Therefore no flux is resolved out of the 
interferometer maps, as expected for these high density 
tracers. 

3. RESULTS 



3.1. The Dense Cloud Morphology 

Continuum subtracted integrated intensity line maps 
of HC3N(5-4), HC3N(10-9) and HC3N(16-15) in IC 342 
are shown in Figure m Fi gure [2] shows the 5-4 transi- 
tion overlaid on the 7 mm continuum image generated 
from offline ch annels. Lo cations of giant molecular cloud 
(CMC) cores ( |Downes et al. 1992; Meier & Turner 2001] 
and the optical clusters (e.g. |Schinnerer et^ al. 2003|\ 
coi npare d with HC3N(5-4) and HCN(l-O) fDowncs et] 
al. 1992), are also shown in Figure I2] Figure I3] shows 
tEIeTTCiN(5-4) and HC3N(16-15) spectra (in flux units) 
taken over the same 2" aperture centered on each cloud. 

HC3N(5-4) emission picks out most clouds seen in 
other dense gas tracers (e.g HCN(l-O). The only la- 
beled cloud not clearly detected in HC3N(5-4) is CMC B, 
the cloud associ ated with the nuclear star-form ing region 
(LiR ~ 10^ Lq; [Becklin et al.|l980| [Turner fc Ho 1983). 
Positions of the GMCs measured in HC3JN are consis- 
tent with those fitted in Ci*0(2-1) to within a beam 
( [Meier fc Turner[[2001[ ). An additional CMC is detected 
in HCgiN (5-4) just south of CMC C3, labeled C4. Unlike 
Ci80(2-1) which peaks at CMC C2, HC3N(5-4) emission 
peaks farther north, towards CMC CI, at a distance of 
105 pc from the nucleus suggesting changes in excita- 
tion across CMC C. GMCs A and D' are resolved into 
two components, but we do not discuss them as sepa- 
rate entities given the lower signal-to-noise, other than 
to state that there is a clear difference in velocity cen- 
troid between the 5-4 and 16-15 transitions (e.g. Figure 
indicating there is likely an excitation gradient across 
TjMC a. No emission > 2a is detected towards the weak 
CO (2-1 ) feature associated wit h the central nuclear star 
cluster ( Schinnerer et al.[|2003 ). 

The J=16-15 line intensity is comparable to that of 
J=5-4 line in all clouds. While dominated by CMC 
C, HC3N(16-15) is detected (or tentatively detected) to- 
wards all clouds, except GMC D. Fluxes for GMCs D and 
D' are quite uncertain since they are just inside the half 
power point of the PdBI primary beam. That GMC C is 
much brighter than the other clouds in both transitions 
indicates large quantities of dense gas are present in this 
cloud. HC3N(16-15) favors C2 over CI. GMC B, while 
still faint in an absolute sense, is significantly brighter 
in HC3N(16-15) relative to 5-4. This is not unexpected 
given the higher excitation requirements of the 16-15 line 
and that both B and C2 have 7 mm continuum sources 
associated with massive star formation. 



3.2. Gas Excitation and the HC3N A J Line Ratios 

Comparisons of the HC3N (5-4), (10-9) and (16-15) 
maps provide a chance to establish the excitation of the 
densest molecular cloud gas component. Line intensi- 
ties from the J=5-4 and J=16-15 lines were measured 
over 2" apertures centered on each of the GMCs. Table 
[l[ records gaussian fits to each spectrum along with in- 
tegrated line intensities. Peak antenna temperature ra- 
tios are calculated for HC3N(16-15)/HC3N(5-4), here- 
after denoted R16/5. R16/5 ranges from less than 0.06 
up to ~0.5, with GMC B having the highest value (Ta- 
ble [l]). The 5-4 and 16-15 transitions generally bracket 
the peak of the level populations, so we achieve good 
constraints on gas excitation. LTE excitation tempera- 
tures, Tex, implied by R16/5 range from <10 K to >18 
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K (where we have neglected Tcmb in this determination) . 
Excitation is lowest towards GMCs D, D' and E. These 
excitation temperatures T^x are similar to those found 
from the presumably much less dense gas traced in C^^O 



(Table 2; |Meier fc Turner|2001 ). Only the CMC C clouds 
have signihcantly higher m HC3N — towards CI - 
C3 Te2:(HC3N) are about a factor of two greater than 

The resolution of the HC3N(10-9) data is significantly 
lower than it is for the 5-4 and 16-15 lines. To compare 
J=5-4 and J=10-9 line intensities, the 5-4 data were con- 
volved to the resolutio n of the (10-9) data (5."9 x 5."1 
Meier & Turner 2001), then integrated intensity ratios, 
hereatter K10/5, were sampled at the locations of Ri6/5- 
Though at lower resolution than Rig/5, we make the ap- 
proximation that Rio/5 does not change on these sub- 
GMC scales. While leading to larger uncertainties, this 
provides a way to include all three transitions in model- 
ing dense gas excitation at very high resolution. i?io/5 
range from 0.25 to 1.4 (Table [l]). GMC B has the high- 
est value of 1.4. The remainder of the GMCs have ratios 
of 0.3 < -Rio/5 < 0.6. Excitation temperatures implied 
by these ratios range from T^x = 6-16 K, consistent 
with those derived from Rie/5 separately. Tex derived 
from Rio/5 are similar to those derived from C^*0. The 
only exception here is GMC A, the clo ud where PDRs 
(Phot on-Dominated Regions) dominate (Meier & Turner 
2005| . Towards GMC A T^xiC^^O) and TexiRie/s) are 
12 - 13 K, twice that of Te^(Rio/5). 

Before modeling the densities implied by the line ratios, 
we test whether IR pumping can be responsible for the 
observed excitation. IR pumping is important if. 



(1) 



where '"^''Bui is the Einstein B^i of the corresponding 
1/^ = 1 vibrational transitions at ~ 45 /im, /iy(45 fim) 
is the 45 /im intensity as seen by the HC3N molecules 
and ^°*Aui is the Einstein Aui of the rotational state 
(e.g. Costagliola & Aalto 2010). It is extremely diffi- 
cult to estimate the applicable 45 iim IR intensity, but 
it is expected to be most intense towards the starburst 
GMC (B). A detailed assessment of IR pumping must 
await high resolution MIR maps, but we constrain /^(45 
lim) in several ways . First, we take the 45 /im fiux from 
Brandl et al. ( 2006 ) and scale it by the fraction of total 



20 /im flux t hat comes from within 2.1 of the starburst 
as found by Becklin et al. (1980) and then average over 
that aperture, b'oi the average (45 /xm) calculated this 
way, '"''Bui -fi/(45 /im) is 10"'' (lO'^-^) times too low 
to pump the 16-15 (5-4) transition. Alternatively if we 
(very conservatively) take the total MIR luminosity from 

the central 30 and assume it comes from a bl ackbody 
of the obse rved color temperature (e.g. ^50 K; [Becklin 
et al. 1980) then 7^(45 /im) is still at least an order of 
magnitude too low to meet the inequality in eq. [l] for 
both transition. IR pumping rates only become compa- 
rable to ^°*Aui for the 5-4 transition if the total IR lu- 
minosity originates from a ~2.5 pc source with a source 
temperature of ZlOO K. We conclude that IR pumping 
is not important for the 16-15 transition of HC3N in any 
reasonable geometry of the IR field. For the 5-4 transi- 
tion to be sensitive to IR pumping, the IR source must 



be warm, opaque and extremely compact. Therefore IR 
pumping is neglected for all clouds in the following dis- 
cussion. 

3.3. Physical Conditions of IC 342's Dense Gas — 
LVG Modeling 

The values of the excitation temperature constrain 
the density and kinetic temperatures, n^/^ and Tj,, and 
the physical conditions of the clouds driving the excita- 
tion. A series of Large Velocity Gradient (LVG) radiative 
transfer models were run to predict the observed intensi- 
ties and line ratios for a given n//^ , T^ and filling f actor. 



fa — ^5 



of the dense component (e.g., 



Van- 



s are 



den Bout et al. 1983 ). Single component LVG mode 
instructive, particularly when lines are optically thin, as 
is the case f o r HC3 N. The LVG model used is that of 
adapted to HC3N, with levels up 



Meier et al. 
to J =20 inc 



& Chapman 



mm, 

u ded. Collision coefficients are from IGreenI 
(|1978|). A range of densities, n//^ — 10^- 



10' 



cm and kinetic temperatures, T^ = 0-100 K, was 
explored. HC3N column densities based on LTE excita- 
tion (Table [1]), are calculated at 2 resolution from: 



Nmol — 



3fcQe^"/''^- 
Sn^vSuHJ-lgK^gi^ 



(2) 



Lafferty 



( |T978 



& Lovas 
.(HC3N) from TSE^ 
found to be, X{HC^N) 



I 



using molecular data of 
IIC3N(5-4) intensities, and' 
IIC3N abundances are 

10-91 _ 10-8-5, with the highest values towards CI 
and D. While uncertain these a bunda nces agree with 
those found in Meier & Turner (2005) an d are typical 
of Galactic center HC3JN abu ndances (e.g. Morris et al 



1976 de Vicente et al. 20001 and good enough tor con 
straining A/du/dr. I'he ratio of cloud linewidth to core 
size is '-^ 1 — 3 km s^^pc~^ for the GMCs. Therefore, we 
adopt a standard model abundance per velocity gradient 
of X/dv/dr = IQ-^ km s"^, but run models for values of 
X/dv/dr 10^" - 10-9. 

Antenna temperatures are sensitive to the unknown 
filling factor. In these extragalactic observations, the 
beam corresponds to scales large compared to cloud 
structure, and hence filling factors are not directly 
known. To first order, the line ratios, -R10/5 and R16/5 are 
independent of filling factor if we assume that HC3N(5- 
4), (10-9) and (16-15) originate in the same gas. So 
model ratios are compared to the observed data to con- 
strain parameter space. For the parameter space implied 
by the line ratios model brightness temperatures are de- 
termined. A comparison of the model brightness tem- 
perature to the observed brightness temperature sets the 
required areal filling factors for that solution. 

3.3.1. Physical Conditions of the Dense Component 

Figure |4] displays the results of the LVG modeling. Ac- 
ceptable ±1(7 [Tk^n-H^) parameter spaces are shown for 
each line ratio (R10/5, thick solid gray line; Rig/s, thick 
solid black line) . Also shown in Figure [4] are the accept- 
able solutions obtained fr om the C^^Pi^-l) IG^^ O(l-Q) 
line ratio (thin gray lines; Meier fc Turner|2001 1. Figure 
[5] and Table [1] display the model fiux versus upper J state 
of the line for the T^ = 10 K, 30 K and 50 K solutions 
(also 70 K solutions for GMC B). 
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HC3N(5-4) has a critical density of 10 cm 

and an upper energy state of 6.55 K. HC3N(16-15) has 
a critical density, ricr ^ 5 x 10^ cm~'^ and an upper en- 
ergy state of 59 K. HC3N(10-9) values are intermediate. 
Therefore i?io/5 and Rie/s are sensitive probes of njy^ 
between the range of 10^ and 10^ cm~'^. Observed ratios 
constrain njy^ to ~ ±0.1 dex at a given T^, and to ~ ±0.3 
dex for all modeled T^ > 20 K. For Tfc > 40 K, ratios 
are largely insensitive to T^ (curves are horizontal). The 
LVG models do not constrain density when T/j is low. 
Thus to narrow the range of possible solutions requires 
two or more line ratios, or an external constraint on one 
of the axes. Kinetic temperatures of the dense compo- 
nent at 2 scales in IC 342 remain largely unconstrained 
to date. One method of constraining the kinetic temper- 
ature is the peak Tf, at this resolution of the optically 
thick C0(2-l) hue. Using this method, Tfc ~ 35-45 K 
towards GM Cs A, B and C and Tfc - 15-20 K for GMC s 
D, D' and E ( [Turner et al. l993[ [Schinnerer et al^|2003[ ). 
These values are not highly discordant from a rcimnute 
resolution Tfc ~ 50 K measurements from NH3 (Ho et al. 



1982 1 or the far-infrared dust temperature o f 42 - 55 K 
( Becklin et al]|1980[ [Rickard fc Harvey||1984| ). However, 



the CO (2-1) traces gas with densities two to three orders 
of magnitude lower than these HC3N transitions, so it is 
not clear that this is the relevant Tfc for the dense cores 
of the GMCs. 

Valid solutions are found for the GMCs, with the ex- 
ception of CMC D', which is the most uncertain due to its 
distance from the center of the field. (No fit is attempted 
for CMC D because both HC3N(5-4) and HC3N(16-15) 
are upper limits or tentative detections.) The solutions 
show modest cloud excitations given the strong nuclear 
star formation. Good agreement is observed between so- 
lutions found using R16/5 and R10/5 when kinetic tem- 
peratures are low (Tfc < 30 K). Given uncertainties in 
line intensities and X{HCzN)/dv/dr, higher tempera- 
tures cannot be ruled out. While Tfc itself is less well 
constrained the combination nH2T^k is well determined. 
Since line brightnesses are the measured quantity, as gas 
temperatures are raised, densities or filling factors must 
decrease to compensate. The nature of the solutions are 
such that the densities decrease more than the filling fac- 
tors (Table |4]). 

GMCs A, CI and C3 are best fit with nj^^ ~ 
104.7-5.2 ^^-3 _ 20 K. For Tfc = 50 K derived 

densities drop to n/f^ ^ 104-3-4.7 gj^^-s^ while they drop 
to UH^ ~ for Tfc = 100 K. Changes in 

X/dv/dr do not strongly influence the derived solutions. 
Towards these clouds filling factors are fa —0.03 - 0.1. 
For comparison, a 1 pc^ cloud would have fa — 0.0012 
for the aperture size and IC 342 's distance. Therefore 
these GMCs appear to have a few dozen dense clumps 
similar to the larger clumps found in Galactic GMCs (e.g. 
Myers fc Benson[1983l [Zinchenko et al.|[1998| ). CMC C2 
appears to have similar densities but slightly elevated 
temperatures, Tfc ~ 30K. In GMCs D and D' densi- 
ties are at least a factor of four lower than the other 
clouds. Little can be said about the kinetic tempera- 
tures of these GMCs. CMC D' is the only location with 
statistically discrepant solutions from R10/5 and Ri6/5- 
Here the (16-15) line is somewhat brighter than would 



be expected from the R10/5. 

The starburst GMC, B, is markedly difii'erent from the 
others. The combination of a low absolute HC3N bright- 
ness temperature and high ratios, -R10/5 and R16/5, re- 
quire higher densities of nH2^^^^ cm~^ and a smaller 
{fa <0.02) filling factor. Both line ratios seem to indi- 
cate a kinetic temperature of more than Tfc ~ 40 K, 
and perhaps significantly higher, averaged over our 30 
pc beam. The dense gas in GMC B, which is closest 
to the active star-forming region and IR source, is more 
compact, denser, and more highly excited than the other 
clouds. 

For our own Galactic Center it has long been argued 
that molecular clouds must be unusually d ense, at least 



19891 
clou 



10 cm to withstand tidal forces (eg. Stark et al 



If one assumes that GMC C, the most well defined 



is a point mass in a spherical mass potential of 
the galaxy, its Roche limit would be at a radius of 25 pc. 
The radius of the combined C1-C2-C3 complex is >25 
pc. In IC 342, the proximity to the nucleus w here the 
strong noncircular velocity field (Turner fc Hurt 1992| 
together wi th feedback from the starburst (Schinnerer 
et al. 2008) suggests that these clumps have densities 



large enough to maintain their identity but likely will 
not remain gravitationally bound to each other. 

3.3.2. Comparisons Between HC3N and C^^O Physical 
Conditions 

For aU GMCs (including GMC B), Te^(HC3N) ~ 
Tex{G^^O). If a model in which all the H2 exists in one 
uniform component, then the similarity of T^x from both 
HC3N and C^^O implies that the densities of the molecu- 
lar clouds are high enough (> 10^'^ cm""^) to thermalize 
both C^^O and HC3N across the nucle us. (LVG solu- 
tions f or C^®0 J=2-l and 1-0 lines from ' Meier fc Turnerl 
(2001 1 are shown with the HC3N solutions in Figure E ) 
In this monolithic model the low observed T^x of 10-2(J~K 
demand that the thermalized clouds must be quite cool. 
These temperatures are similar to those of dark clouds 
in the disk of our Galaxy, which is somewhat surprising 
given the elevated star formation rate in the nucleus of 
IC 342. 

However this model of a monolithic, dense, cold ISM 
seen in both HC3N and C^®0 runs into problems with 
total mass constraints. If the ISM is uniformly this 
dense and cold then the total dense gas cloud mass im- 
plied by the required densities and filling factors become 
very large, greater than permitted based on the optically 
thin C^^O line emission. In Table E] the total mass of 
dense gas, Mden, is approximated from the LVG solu- 
tions assuming Mrfe„ rTi//^ ri// ^ -R'^ , with R defined as 
R — V faAbeam- Table [l] lists the total mole cular gas 
mass, , e stimated from C^®0(2-1) data (see Meier fc 
Turner|2001 ) assuming Mj^^ = mH2NH2^beam- One can 
see that Mden is typically larger than Mh2 for Tfc < 30 
K and diverges rapidly as Tfc drops towards the LVG fa- 
vored values. In short if GMCs have uniformly such high 
densities and cold temperatures then the clouds would 
contain too much H2 mass for what is observed in C^^O. 
To match the line ratios while not violating mass con- 
straints requires a multi-component dense ISM. A possi- 
ble two component model of the dense gas is discussed 
in SlOl 
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4. DENSE GAS AND THE CO CONVERSION FACTOR 

A consideration of the conversion factor between CO 
intensity and H2 column density, Xco^ in the nuclear 
region of IC 342 suggests that the model of a single, 
high density and relatively cool ISM is not consistent 
with observations of CO isotopologues. Moreover these 
clouds are unlikely to resemble Galactic disk GMCs in 
their internal structure and dynamics. 

4.1. The Single Component Model and Xco 

The well known Galactic relation between CO inten- 
sity and H2 column density, Xco — Ico/X{H2), can 
be explained if GMCs consist of optically thick (iii CO) 

La rson l lQSl 



turbulent clumps in virial equilibrium (e. 



Solomon et al. 1987 Scoville & Sanders 111957). b'or the 



clumps to emit m MCsiN densities must be large enough 
so that Tf, can legitimately be approximated by T^. If 
we adopt this model, then Xco is approximated by: 



4.2. A Possible Two Component Model and Xco 

The lack of an apparent connection between V'^^n/Tf^ 
of the dense component and the observed conversion fac- 
tor indicates that the clouds in the nucleus of IC 342 
cannot be treated as simple virialized collections of uni- 
form density clumps. One might expect that a large frac- 
tion of the C^^O emission could originate from a moder- 
ate density component that is distinct from the denser, 
HCsN-emitting gas. (Note that ^■^CO traces a distinct 
component more d iffuse than that traced both by C^^O 
( |Meier et al.|[2000l ) and HC3N.) 

As a tirst order extension to the basic virialized clumps 
model, we imagine two spatially well-mixed sets of 
clumps; one low density C^^O emitting, ni, and one high 
density, n^, that emits both C^^O and HC3N. Assuming 
the fraction, by number, of clumps with high density is 
r, then if both clumps have the same T^ the virialized 
clumps Xco relation becomes: 



Xco - 0.84 X 10^" 



fa ^'Tfc 



cm-2(Kkms-i)-i (3) ^^^^^Xco ^ 0.84 x 10^ 



y/4> n 

fa Tk 



cm ^ (K km s 



where 4> is the volume filling factor and "^'n is the densit y 
of the clumps (Sakamoto 1996 Maloney & Black 1988). 
From the bri ghtness of the CO iso topic lines, ^'^/g ~ 1 
on 2" scales ( Meier fc Turner|2001 ). The conversion fac- 
tor Xco within the central tew hundred pc of IC 342 
has been determined to be 0.6 x 10^° cm^^ (K kms^^)" 



(Meier & Turner 2001). If the monolithic model pre- 
sented in the previous section is correct and clouds are 
virialized clumps in equilibrium, then V'^'n/Tfe must be 
^0.75 to match the observed conversion factor. 
Table [1] records \A%i/Tfc for each of the LVG solutions 

across the nucleus of IC 342. For all solutions, "v/^/Tfc 
is greater than the required value, especially at low Tfc. 
For the single component LVG favored Tj. of ~10 - 20 
K towards the more quiescent clouds, V'^'n/Tfe ~5-75, 
implying ^'^^'^^Xco should be several times larger than 
*^^Xco- However Xco in the nucleus of IC 342 is 3-4 
times lower than the Galactic conversion factor, based 
on both optic ally thin isoto pologues of CO as well as 
dust emission (Meier & Turner 2001). Lower conversion 
factors appear to be the norm for the nuclei o f gas-rich 
st ar form ing galaxies, including our own (eg. [Smith et 
ar]|l99T| ISc oville et al. 1997; Dahme n et al.||l9fe| JHarri- 
son et al. |1999j|Meier fc Turner|2004| [Meier et alj|2008p . 



■ ■■ ^1- 

In the nuclear region ot IC 342 dense gas kinetic tem- 
peratures do not appear to be higher than the Galactic 
disk value by an amount large enough to offset the ob- 
served increase in density over typical Galactic dis k-like 
clouds ; Tfe >100 K would be required. Recently |Wall| 
(2007) completed a more generalized treatment of the 



conversion factor including radiative transfer and con- 
cludes that the true exponential dependences of n^/j and 
Tfc are weaker than 0.5 and -1, respectively. However in 
all cases modelled, the dependences remain positive for 
Uh^ and negative for Tfc. 

In summary, it is not possible to reconcile bright HC3N 
emission from uniformly cool, dense gas (low Tfc, high 
and high fa) with the known total amount of H2 
present if virialized constant density clumps are adopted. 



(4) 

where n = (1 — r)rt; + Tn^. The ratio of dense gas to 
total H2 becomes, Mden/'^tot ^ rn,i/((l - T)ni + Tnu) 
— Tufi/n. For a given Tfc, ni is chosen to match 
the Ci^O LVG solution ( [Meier fc Turner||200l| ), while 
nil is chosen to match the HCaiN LVG solutions, and 

r = (f>HC3N /(t>co — fa^^- Table [l] displays the adopted 
C^^O LVG solutions, along with the new \/n/Tk and 
Mden/Mtot for the two component model, assuming, for 
simplicity, that Tfc is the same for both components. 
Vn/Tfc and hence Xco are decreased from the simple 
one component model by factors of at least four. This 
simple extension results in closer agreement between the 
observed and calculated values of Xco (especially when 
including the somewhat super- virial line widths), while 
maintaining most of the mass in the dense (low filling 
factor) component. 

It is almost certainly the case that this extension is 
an oversimplification. In reality, we expect the clouds to 
exhibit a continuum of densities. However the LVG mod- 
eling demonstrates that at least these three components 
are required to match the multi-transition observations 
of i^CO, i^CO, C^^O and HC3N. This is consistent with 
conclusions from rece nt single-dish modeli ng of higher J 
transitions of HC3N (Aladro et al. 2011), but we find 
multiple components are required to match intensities 
not only between GMCs but within individual GMCs. 

5. HC3N EMITTING DENSE GAS AND STAR FORMATION 
5.1. HC3N versus HCN(1~0) 
HCN(l-O) emission is the workhorse f or relating quan- 



tities of dense gas to star formation (e.g. Gao fc Solomon 
2004). It is interesting to compare conclusions about 
the dense component from HC3N with those of the more 
commonly used HCN(l-O). H CN(l-O) has b een mapped 
at similar spatial resolution by Downes et al. ( 1992). Fig- 
ure [2] compares HC3N(5-4) to HCiN(l-O). While HC3N 
generally traces the same dense GMCs seen in HCN their 
relative brightnesses are rather different. In HCN (1-0) 
CMC A, B and C are all within 10 % of the same bright- 
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ness and GMCs D and E are 100 % and 50 % weaker, re- 
spectively. Whereas in HC3N(5-4) and (16-15), C dom- 
inates, B is nearly absent and A is not significantly dif- 
ferent from D and E. Comparisons with HC3N clearly 
demonstrate that there is larger variations in dense gas 
properties than HCN indicates. The dominate difference 
is enhanced HCN towards the starburst and GMC A rela- 
tive to GMC C. Unlike the HC3N transitions, HCN(l-O) 
is optic ally thick and has sli ghtly larger (25 - 50 %) fiUing 
factors ( Downes et al.|1992 |. As kinetic temperatures in- 
crease optically thick transitions brighten more rapidly 
than (lower excitation) optically thin transitions. Hence 
it is expected that the HCN emission should favor some- 
what warmer dense gas. Likely this effect results in the 
much brighter relative HCN(l-O) intensities towards the 
starburst. The relative enhancement towards GMC A is 
less clear. Ho wever, this clo ud is known to be dominated 
by PDRs ( jMeier fc Turner p305) , strongly in fluenced by 
mechanical feedback from the nuclear cluster ( Schinnerer 



et al. 2008 1 and has a complic ated HC3N temperature 
and velocity structure (Sections 3.1 and 3.2). This cloud 
must have a complex density and temperature structure, 
potentially with a warmer intermediate density medium 
(Table IT]) . Moreover, HCN abundances can be elevated 
in PDRs. 

We conclude that while HCN(l-O) does a good job 
locating the dense gas but it does a poorer job tracing 
small scale variation in the properties of the dense gas 
when a mix of strong star formation and quiescent gas are 
present. It is expected that such effects should become 
more important as specific star formation rates increase 
and spatial scales decrease. 

5.2. HC3N versus Star Formation 

The brightest HC3N emission and the brightest star- 
forming regions do not coincide. Most of the current star 
formation, traced by bright infrared and thermal radio 



conti nuum emission ( [Becklin et al.| 1980 Turner & Ho 
1983), is situated about 50 pc to tne southwest of the 
dynamical center, in the vicinity of GMC B. The bright- 
est HC3N emission by far is on the northeast side of the 
nucleus, in the northern molecular arm centered at GMC 
C2. This region has free-free emission amounting to only 
a third of the brightness of the strongest radio source, as- 
sociated with GMC B. The faintness, in absolute terms, 
of HC3N(16-15) towards GMC B is unexpected and im- 
portant reflecting a much lower areal filling of warm, very 
dense gas here. Either the number of cloud clumps or 
their size is small relative to the other GMCs. 

On the other hand, the excitation of HC3N does reflect 
the presence of young forming stars. HC3N(16-15) and 
HC3N(10-9) are relatively brighter towards the current 
starburst. This suggests that the absence of a correla- 
tion between HC3N(5-4) intensity and star formation is 
partly due to depopulation of the lower energy transi- 
tions. GMCs that are faint in the HC3N J=16-15 tran- 
sition and not associated with strong star formation show 
up well at the lower J transitions of HC3N, as would be 
expected for gas of lower excitation. Clearly there is a 
large amount of dense gas currently not actively forming 
stars, that shows up in the low excitation transitions of 
HC3N. This is consistent with the fact that both HNC(1- 
0) and N2H+(l-0), generally considered dense quiescent 
gas tracers, are found to be very bright towards IC 342's 



nucleus ( [Meier fc Turner|[2005| ) . HC3N(5-4) appears to 
be an excellent extragalactic probe of the dense, quies- 
cent molecular gas component not yet involved in the 
current starburst. 

To quantitatively compare star formation with dense 
gas, a star formation rate is derived from the 45 GHz 
continuum flux (spectral indices measured between 2 cm 
and 7 mm demonstrate that the vast majority of flux to- 
ward GMCs B and C is thermal Bremsstrahlung) . The 
star formation rate is then compared with Mden (Ta- 
ble [1]) to estimate a dense gas star formation efficiency, 
SFEden = SFR{Mq yr-i)/A'W(M0). Dense gas de- 
pletion timescales, Tden= l/SFEde„, are also computed. 
In Table [T] the ratio of the observed Lyman continuum 



ioniza tion rate (e.g. Meier & Turner 2001 Tsai et al 
2006 corrected for distance) to M^en derived from the 
LVG analysis is reported 

GMCs NLyc/Mden ^ 1-6 - 4.0 X 10'"S 

wards GMC B this ratio is 10 - 30 times larger. If one 
adopts the conversion fro m N^^^ to st ar formation rate of 

Kennicutt|[r998i ), then SEE 



Towards the non-starburst 
45=-iM0\ while to- 



lO"^'^ Mq yr'^ s- 

for the starburst is '-^ 4 x 10~' yr~^, or dense gas de- 
pletion times of r^en — 2-3 Myr! This highly enhanced 
SEEden is a direct consequence of the faint HC3N emis- 
sion here. Even towards the non-starburst clouds SEE^en 
are 1.6— 4.0x 10~* yr~^. These efficiencies are sufficiently 
short that they imply dense gas consumption timescales 
that are non-negligible fractions of the expected GMC 
lifetimes. 

The dense SEE is rather high across the nucleus, but 
the extreme valu e towards GMC B is remarkable. Meier] 
& Turner (2001) argued that intense star formation is 
suppressed along the spiral arms being triggered when 
the inflowing molecular gas collides with the inner ring 
molecular gas. Therefore it is reasonable that in a rel- 
ative sense SEE^en is lower away from the central ring. 
However this leaves unexplained why GMC B's SEE^en 
is so much larger than GMC C, though its positions at 
the arm / inner ring intersection is the same. We sug- 
gest this is a sign of the evolution of star formation across 
the nucleus that is impacted by radiative and mechanical 
feed-back from within the molecular cloud. 



5.3. Destruction 



' Dispersal of Dense Gas with 
Starburst Age 



A possible cause of the different SEE^en between GMC 
C and B is that we are observing the clouds at high 
enough spatial resolution to begin to identify the chang- 
ing internal structure of the clouds in the presence of the 
starburst. Over the lifetime of a cloud SFEs vary. In the 
earliest stages of a star formation episode SEEs will ap- 
pear low because elevated star formation rates have yet 
to convert the bulk of the molecular material to stars. 
Towards the final stages of a GMCs evolution instanta- 
neous SEEs appear to increase dramatically as the cloud 
clumps are consumed, destroyed or dispersed. So ob- 
served instantaneous SEEs are expected to vary widely 
throughout the lifetime of an individual GMC and rel- 
ative to lifetime averaged SEEs typically considered in 
extragalactic studies. 

If starburst B is a few Myrs more evolved than the 
other GMCs, especially the dynamical similar C, then 
we may be witnessing the consumption, dispersal or de- 
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struction phase of the remaining dense clumps in the 
presence of the expanding HII region. The magnitude of 
the dense gas consumption times for B are indeed shorter 
than the lifetime of the GMC. Several lines of evidence 
suggest that the (weaker) star formation towards C2 may 
be at a somewhat earlier phase. These include less ex- 
tended HII regions (e.g. Tsai et al. 2 006 ) , brig ht hot core 
like speci es CH3OH (jMeier fc Turner 
jMonter o-Castafio et al.||200tip and CU _ 
et al. | 199ll, and more mm dust continuum emission (e.g. 
Meier fc 'iuT^|200T| ). 

in this context it is interesting to compare the ther- 
mal pressure of the starburst HII region with that of the 
dense clumps along the same line of sight. Assuming a 
Stromgren sphere of R~3 pc, N^yc = 1.1 x 10^^ s~^ and 
Te ~ 8000 K, parameters determined for the main star- 
burst HII ( Tsai et al.|2006 ) , the thermal pressure of such 



t)-5) (e.g. Harris 



an HII region would be rigTe = 10 cm"'' K. This value 
is equal given the uncertainties to na-iYK — lO^-^''"^ ''^ 
cm~^ K for the dense component of GMC B, hinting at 
pressure balance between the dense clumps and the HII 
region. In addition, the filling factor of the HII region 
is larger than implied by the HC3N LVG analysis. So it 
is possible that HC3N emission towards GMC B comes 
from dense clumps embedded within the HII region. 

The above analysis suggests the following physical pic- 
ture for the faint HC3N emission towards the starburst. 
The starburst towards B is more evolved. The HII re- 
gion at GMC B has had time to expand, destroying or 
dispersing the dense gas, which is now in the form of 
smaller clumps and/or more diffuse gas. Clumps that 
remain there must have high pressure to survive. Near 
the younger star forming cloud GMC C (particularly to- 
ward GMC C2), the HII regions may just be developing, 
and have not had time to disperse the clouds. The dense 
clumps here would be more abundant and still present a 
hot core-like chemistry. The clouds further from the cen- 
tral are on average less dense and at the current epoch 
remain largely quiescent, except possibly D', where the 
second high excitation com ponent could be assoc iated 
with the presence of shocks ( [Meier fc Turner||2005[ ) . 



6. CONCLUSIONS 

We have imaged the HC3N J=5-4 line in the nucleus 
of IC 342 with the VLA and the HC3N(16-15) line with 
the PdBI at <2" resolution. These are the first maps 
of these transitions in an external galaxy. We have de- 
tected emission extended along the nuclear "mini-spiral" 
in (5-4) and more concentrated emission in (16-15), with 
relative abundance of X^HC^N) - 1 x lO^^. HC^N 
emission is not tightly correlated with star formation 
strength. Dense gas excitation however, follows star for- 
mation more closely. GMC B, which is weak in all the 
HC3N lines, is relatively stronger in the higher J lines. 

LVG modeling indicates that the HCsN-emitting gas 
has densities of lO'*^^ cm""^. In IC 342, physical condi- 
tions of the densest component are fairly constant away 
from the immediate environment of the starburst, though 
beyond the central ring densities begin to fall. Co i npari- 
son with the C^^O observations of [Meier fc Turner[ ( [2001 1 
reveal excitation temperatures similar to C^^O values in- 



dicating either that the molecular gas is dense and cool 
(Tfc < 20 K) or that there are multiple gas compo- 
nents where the densities and kinetic temperatures of 
each component conspire to give similar overall excita- 
tion. The strong over-prediction of the amount of gas 
mass present if densities are large and temperatures cool, 
favors a multi-component ISM with at least two compo- 
nents beside the diffuse one seen in ^^C0(2-l). The ac- 
tual ISM is likely a continuum of cloud densities with dif- 
ferent densities dominating different tracers. HC3N also 
differs in morphology from HCN(l-O), with HCN(l-O) 
being much brighter towards the starburst. This is fur- 
ther evidence that there are multiple dense gas compo- 
nents. 

Of particular interest, the starburst site (GMC B) ex- 
hibits the largest difference in intensity between HC3N 
(both transitions) and HCN(l-O). The faintness of the 
HC3N here suggests that the brightness of HCN(l-O) is 
not due solely to large quantities of dense gas. A compar- 
ison between the GMCs with the largest star formation 
rates and similar dynamical environments, B and C, hint 
at an explanation. While GMC B shows higher excita- 
tion, the low brightness of this cloud indicates that it is 
composed of a relatively small amount of warm, dense 
clumps. The smaller amount of dense gas at the site of 
the strongest young star formation indicate high star for- 
mation efficiencies in the dense gas. Towards GMC C, 
HC3N, CH3OH and NH3 are more intense and the frac- 
tion of millimeter continuum from dust is higher. This 
indicates that GMC C is in an early, less evolved (hot 
core-like) state. The extreme dense gas star formation 
efficiency observed towards GMC B reflects the fact the 
main burst is in a more evolved state. The dense clumps 
towards the starburst are being dispersed or destroyed 
in the presence of the HII region. The little dense gas 
remaining appears to be in pressure equilibrium with the 
HII region. The larger opacity of HCN(l-O) relative to 
HC3N elevates its brightness temperature in this warm 
gas and lowers its critical density permitting it to remain 
excited in the somewhat lower density component. 

We conclude that EVLA observations of HC3N(5-4) 
can be a powerful probe of dense, quiescent molecular 
gas in galaxies, and when combined with high resolution 
imaging of the higher J transitions of HC3N with current 
and upcoming mm interferometers (like ALMA) provide 
tight constraints on dense molecular gas properties in 
stronger or more widespread starbursts, where changes 
like those localized to GMC B are expected to permeate 
much of the ISM. 
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TABLE 1 
HC3N Measurements 



GMC I(5-4)=>- Spfc(5-4) Vo{5-4) At;(5-4) 1(16-15)=' Spfc(16-15) Do(16-15) Av{16-15) 





(K km s-i) 


(mJy) 


(km s-i) 


(km s-i) 


(K km s-i) 


(mJy) 


(km s-1) 


(km s-i 


A 


13±2.6 


5.0±0.7 


22±2 


29±5 


2.4±0.5 


7.0±1.3 


16±3 


44±10 


B 


<5.3 


<1.0 






2.4±0.5 


5.2±1.5 


19±3 


28±10 


01 


44±2.6 


6.6±0.5 


48±1 


39±3 


8.5±0.5 


13±1.1 


47±1 


43±4 


02 


28±2.6 
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10±0.5 
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45±1 
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03 
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45±1 
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D 
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1.5±0.5 
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<1.0 


<2.5 






D' 


20±2.6 


3.8±0.5 


52±2 


40±7 


<1.0 


2.3±1.1 


57±8 


44±26 


E 


14±2.6 


4.4±1 


12±2 


28±7 


<1.0 


<2.5 







Note. — Based on spectra from 2 apertures centered on each cloud, except where noted. Uncertainties in the 
temperatures and intensities are the larger of the rms or 10 % absolute calibration uncertainties. Uncertainties 
for the line centroids and widths are the la uncertainties in the gaussian fits. 

" Based on the full resolution data. 
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TABLE 2 
Line Ratios and Excitation 



GMC 




10/5T a 
-1- ex 


-^16/5*' 




ClSOrn c 

-L ex 


N(H2)<= 


X(HC3N) 






(K) 


(K) 


(K) 


(cm^^) 




(Me) 


A 


0.55±0.2 


8.8±1.5 


0.14±0.03 


12±0.6 


13±4 


4(22) 


l(-9) 


4.4(5) 


B 


2.5±1 


37±19 


^0.52 


^18 


19±8 


4(22) 


>8(-10) 


7.6(5) 


CI 


O.SOitO.l 


llitl.O 


0.20±0.03 


13±0.6 


8±3'^ 


7(22) 


3(-9) 


6.8(5) 


C2 


l.libO.l 


14±1.0 


0.33±0.03 


15±0.6 


8±3 


1(23) 


2(-9) 


1.0(6) 


C3 


l.OitO.l 


13±1.3 


0.42±0.07 


17±1.3 


8±3'^ 


5(22) 


l(-9) 


4.9(5) 


D 


;So.i5 


~5 


<0.17 


<13 


6±4 


2(22) 


3(-9) 


6.4(5) 


D' 


~0.24 


~6 


0.061±0.04 


lOibl.8 


~6 


6(22) 


l(-9) 


5.8(5) 


E 


l.litO.3 


14±3.1 


<0.057 


<10 


7±3 


6(22) 


l(-9) 


1.1(6) 



Note. — Based on spectra from 2 apertures centered on each cloud, except where noted. Uncertainties 
in the temperatures and intensities are the larger of the rms or 10 % absolute calibration uncertainties. 
Uncertainties for the line centroids and widths are the Icr uncertainties in the gaussian fits. 

Based on 6" data. 
^ Base d on the full resolution data. 
" From |Meier fc Turner] l |2001| |. 

Assumed to be constant across GMC C. 



TABLE 3 

One Component HC3N LVG Solutions {X/dv/dr = 10"^) 



GMC Tfc 


log(nH2 ) 


fa 




log(nH2Tfc) 


(K) 


(log(cm^'^)) 




(Mo) 


(\og{cm^^ K)) 



A 


10 


5.15 


0.091 


38 


4.7(6)*^ 


11 


6.16 




30 


4.63 


0.051 


6.9 


6.0(5) 


1.4 


6.11 




50 


4.39 


0.067 


3.1 


5.2(5) 


1.2 


6.09 


B 


10 


5.75 


0.015 


75 


1.3(6) 


1.7 


6.75 




30 


5.22 


4.7(-3) 


14 


6.5(4) 


0.084 


6.70 




50 


4.87 


6.0(-3) 


5.5 


4.2(4) 


0.055 


6.57 




70 


4.81 


5.4(-3) 


3.6 


3.1(4) 


0.041 


6.66 


CI 


10 


5.23 


0.13 


41 


9.7(6) 


14 


6.23 




30 


4.69 


0.066 


7.4 


1.0(6) 


1.5 


6.17 




50 


4.50 


0.075 


3.6 


8.0(5) 


1.2 


6.19 


C2 


10 


5.36 


0.12 


48 


1.3(7) 


13 


6.36 




30 


4.80 


0.056 


8.4 


1.0(6) 


1.0 


6.28 




50 


4.64 


0.057 


4.2 


7.3(5) 


0.73 


6.34 


C3 


10 


5.44 


0.076 


53 


7.1(6) 


14 


6.44 




30 


4.86 


0.033 


9.0 


5.3(5) 


1.1 


6.34 




50 


4.70 


0.034 


4.5 


3.8(5) 


0.0.76 


6.40 


D 


10 
30 
50 


<4.49 
<4.10 
<3.88 


>0.040 
>0.040 
>0.060 






<5.49 
<5.58 
<5.58 


D' 


10 
30 
50 














E 


10 


4.90 


0.12 


28 


4.0(6) 


3.6 


5.90 




30 


4.36 


0.096 


5.1 


8.3(5) 


0.75 


5.84 




50 













Mjjense is an estimate of the dense H2 mass implied by those LVG solutions. 
Mden = mjj^njjji?!, with i?e defined as i?^ = faAaper- 
^ Numbers of the form a(b) are equal to a X 10**. 
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TABLE 4 

Two Component HC3N LVG Solutions {X/dv/dr = IQ-^) 



CMC '"'"^Tfe \og{ni) V login) 

(K) {log(cm-3)) (log(cm-3)) 



A 

/\ 


oU 


z.D 


U.Ulo 


o.Ui 


1 1 
i . 1 


U.Dz 


4.0(45) 




50 


2.4 


0.017 


2.82 


0.55 


0.63 


4.0(45 ) 


B 


30 


2.8 


3.2(-4) 


2.83 


0.87 


0.078 


5.2(46) 




50 


2.5 


4.6(-4) 


2.54 


0.37 


0.097 


4.2(46) 




70 


2.4 


4.0(-4) 


2.44 


0.24 


0.092 


4.4(46) 


CI 


30 


2.4 


0.017 


3.03 


1.1 


0.77 


1.9(45) 




50 


2.1 


0.021 


2.90 


0.56 


0.84 


1.8(45) 


C2 


30 


2.4 


0.021 


3.03 


1.1 


0.77 


2.1(45) 




50 


2.1 


0.013 


2.87 


0.54 


0.83 


1.9(45) 


C3 


30 


2.3 


5.0(-3) 


2.80 


0.84 


0.69 


4.1(45) 




50 


2.1 


6.3(-3) 


2.64 


0.42 


0.72 


4.0(45) 


D 


30 
50 


<2.5 
<2.3 












D' 


30 
50 














E 


30 


2.3 


0.030 


2.94 


0.99 


0.78 


1.6(45) 




50 


2.0 











Note. - — Th e low density component comes from the LVG solutions to C^ O | jMeier| 
|fc Turner|2001^ . For sake of simplicity it is assumed that the kinetic temperature ot 
both the low and high density components are the same. See section [4.2| for discussion 
of r and h. 
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Fig. 1. — The continuum subtracted HC3N integrated intensity maps for IC 342. a) The HC3 N(5-4) contoured in ste ps of 7.3 K km s ^ 

1 2005 I contoured in steps 



Meier &c Turner 



(Scr) for a resolution of 1. 95 X 1. 50. h) a continuum subtracted version of HC3N(10— 9) from 
of 1.25 K km s^^ for a beam size of 5. 9 X 5. 1. c) HC3N(16-15) contoured in 3cr steps of 5.1 K km s~^ for a resolution of 1. 83 X 1. 55 
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, i HC3N(5-4) 



O « 



6? 




« 
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J L 
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Fig. 2. — a) The HC3N(5-4) integrated intensity map (black line) for IC 342 overlaid on the 7 mm continuum emission (grayscale) 
extracted from the line free portion of the HC3N(5— 4) datacube. The continuum grayscale ranges from 0.22 mjy bm~^ - 2.2 mjy bm~^. 
HC3N co ntours are as in Figu re [T] b) The HC3(5-4) integrated intensity map overlaid on the HCN(1— 0) emission (grayscale and dashed 
contours; [Downes et al.|l992| . HC3N contours are as in a). The HCN(l-O) gr ayscale runs from 0.50 Jy beam km s~^ to 4 Jy beam km 
s~^ and contours are in steps of 0.75 Jy beam km s^^. Positions of GMC cores ([Meier &: Turner|200l| and optical clusters (e.g. |Schinnerer| 
et al.|2o"o3 with absolute optical positions having absolute uncertainties of ~ 1 ) are labelled. 
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Fig. 3. — The HC3N(16-15) (thin line) and HC3N(5-4) (thick line) spectra. Spectra are summed over 2" apertures centered on each 
position at the locations of the main GMCs. Note that the bright spectral feature at ~-60 km s~^ is H2CO(2o2 — loi)- 
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Fig. 4.— a) The HC3N LVG models of the GMCs in IC 342. The thick black lines delineate the observed ztlfr range of the HC3N(16- 
15)/HC3N(5— 4) line ratio (Ri6/5)i sampled at 2" resolution. The thick gray lines represents the observed ztlo" range of the HC3N(10- 
9)/HC3N(5— 4) line ratio (Rio/s): sampled at 6" resolution. Xi^HC^N) / dv / dr = 10~^ is assumed for the displayed models. The thin 

black lines displays the ztlcr range implied by LVG models of C^**0 for comparison |Meier fc Turner|200l| . Downward (upward) tickmarks 
indicate upper (lower) limits. 







Fig. 5. — HC3N level population diagrams. Fluxes are taken from a 2 box centered on the emission peak. HC3N(10-9) fluxes are 
extrapolated from the (5-4) data using the line ratios in Table [l] For each GMC, fitted LVG solutions are displayed for Tk = 10 (solid 
lines), 30 (dashed lines) and 50 K (dotted lines), except GMC B which includes T/^ = 70 K dot-dashed lines) (TableQ. Only models with 
X{YlC-j,'n)/dv/dr = 10''^ km"! s are shown. 



